format 8 
markov responseblue
*(muT, WF, xFO, muF, CF, laF, xFC, CP, muP, laP, xCP, muC, laC, WP, xCF, xCD, muD, laD) 

func WP(laP,muP) \
1-((laP/muP)^2/(1-(laP/muP)))*(1+(laP/muP)+0.5*((laP/muP)^2)*((2*muP)/(2*muP-laP)))^-1
func WF \
1 - (4*(laF/muF)^4/(24*(4-(laF/muF))))*(1+(laF/muF)+0.5*((laF/muF)^2)+(1/6)*((laF/muF)^3)+(1/24)*((laF/muF)^4)*((4*muF)/(4*muF-laF)))^-1


T F1 muT
F1 Done WF(laF,muF)*xFO*muF
F1 F2 muF*(1-WF(laF,muF))
F1 C xFC*muF*WF(laF,muF)
C P1 xCP*(muC-laC)
C F1 xCF*(muC-laC)
C D xCD*(muC-laC)
P1 C WP(laP,muP)*muP
P1 P2 muP*(1-WP(laP,muP))
D C (muD-laD)
F2 Done xFO*(CF*muF-laF)
F2 C xFC*(CF*muF-laF)
P2 C (CP*muP-laP)
end


end

bind 

muT 0.2
muF .375
muC 1
muD .2
muP .05

WF 0.3
WP  0.4
CF 4
CP 2

xFO 0.5
xFC 0.5
xCF 0.46
xCD 0.33
xCP 0.21

laF 1/15
laC (0.5)/(0.46*15)
laD (0.5*0.33)/(.46*15)
laP (0.5*0.21)/(0.46*15)

end

expr mean(responseblue)
expr variance(responseblue)

cdf(responseblue)


